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S^ Very high energy (~ TeV) 7-rays from blazars are attenuated by photons from 

the extragalactic background light (EBL). Observations of blazars can therefore 
provide an ideal opportunity for determining the EBL intensity if their intrinsic 
spectrum is known. Conversely, knowledge of the EBL intensity can be used to 
determine the intrinsic blazar spectrum. Unfortunately, neither the EBL inten- 
sity nor the intrinsic blazar spectrum is known with high enough precision to 
accurately derive one quantity from the other. In this paper we use the most 
recent data on the EBL to construct twelve different realizations representing 
all possible permutations between EBL limits and the detections in the different 
wavelength regions. These realizations explore a significantly larger range of al- 
lowable EBL spectra than any previous studies. We show that these realizations 
can be used to explore the effects of the EBL on the inferred spectra of blazars. 
Concentrating on the two relatively nearby (z w 0.03) blazars Mrk 421 and 501, 
we derive their intrinsic spectra and peak 7-ray energies for the different EBL 
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realizations. Some EBL spectra give rise to "unphysical" intrinsic blazar spectra, 
characterized by an exponential rise at high TeV energies. We use the F-test to 
show that some of these exponential rises are statistically insignificant. Conse- 
quently, statements regarding the existence of a "IR background- TeV 7-ray crisis" 
are unfounded on the basis of our present knowledge of the EBL. EBL spectra 
that do give rise to unphysical blazar spectra are regarded as invalid realizations 
of the EBL. Those that do not, define new constraints on the EBL spectrum, 
and are used to derive new limits on the intensity and the peak 7-ray energy of 
these two blazars. In particular, we derive an upper limit of ~ 15 nW m~ 2 sr _1 
on the 60 /iin EBL intensity, and find the peak energies of the two Mrk 421 and 
Mrk 501 sources to be very similar, between 0.5-1.2 TeV and 0.8-2.5 TeV, re- 
spectively. We also show that the intrinsic spectrum of Mrk 421 during a period 
of intense flaring activity has a peak energy that seems to shift to higher energies 
at higher flux states. Finally, we show that most EBL realizations show that 
the spectrum of the most distant blazar H1426+428 (z = 0.129) peaks between 
~ 1 and 5 TeV, with some suggesting peaks below 400 GeV and above 10 TeV. 
These results provide important constraints on the different particle acceleration 
mechanisms and the generation of the 7-ray emission. Finally we also show that 
uncertainties in the absolute calibration of the 7-ray energies with atmospheric 
Cherenkov telescopes have an important impact on the intrinsic blazar spectra. 
These systematic uncertainties will be improved with the next generations of 
telescopes which will also cover a wider range of 7-ray energies, providing further 
insights into the intrinsic spectrum of TeV blazars. 

Subject headings: BL Lacertae objects: individual (Markarian 421, Markarian 
501, H1426+428) - galaxies: active - gamma rays: observations - cosmology: 
diffuse radiation - infrared: general 



1. INTRODUCTION 

Observations of blazars at TeV 7-ray energies provide a view of the high energy tail of 
blazar spectra, and are instrumental for probing non-thermal phenomena in jets of active 
galactic nuclei. These spectra provide important constraints on the particle acceleration 
mechanism that generates 7-ray photons in the jets of active galactic nuclei, the cooling 
processes by which particles dissipate their energy into electromagnetic radiation, and on 
the astrophysical conditions in the acceleration region. TeV observations therefore might 
offer a powerful way for discriminating between different 7-ray emission models, recently 
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reviewed by Coppi (2003). 

To the dismay of 7-ray astronomers, high energy 7-rays traveling cosmological distances 
are attenuated en route to Earth by 7 + 7 — > e + + e~ interactions with photons from the 
extragalactic background light (EBL; Gould & Schreder 1967; Stecker, de Jager, & Salamon 
1992). Consequently, particle acceleration mechanisms and 7-ray production models can 
only be tested after the observed blazar is corrected for this attenuation. 

However, to the delight of infrared (IR) astronomers, the attenuation of TeV 7-rays can 
be used to constrain the intensity and spectrum of the EBL in the wavelength region where 
its direct measurement is most difficult. Defined here as the diffuse background light in the 
0.1 to 1000 /im wavelength regime, the EBL is the second most important radiative energy 
density in the universe after the cosmic microwave background. It is the repository of all 
radiative energy inputs in the universe since the epoch of decoupling (Partridge & Peebles 
1967; Hauser & Dwek 2001). It therefore contains important information on the release of 
nuclear, gravitational, and possible exotic forms of energy throughout most of the history of 
the universe. Unfortunately, direct measurement of the EBL in the ~ 5 to 60 /im wavelength 
region is greatly hampered by foreground emission from interplanetary dust (Hauser et al. 
1998; Kelsall et al. 1998). Consequently, TeV 7— ray observations may be the only way to 
constrain or determine the EBL in this wavelength regime. 

Early studies (e.g. Stecker, de Jager, & Salamon 1992; Dwek & Slavin 1994; Biller et 
al. 1995) assumed that the intrinsic spectrum of blazars is a power law in energy E of the 
form dN/dE oc E~ a , and searched for a "break" in the spectrum that could be attributed to 
attenuation by the EBL. In a more recent paper, Costamante, Aharonian, Horns, & Ghisellini 
(2004) used assumptions that reach beyond these early attempts relying on a "break" in the 
spectrum: they require that the spectrum be "smooth", to derive limits on the EBL. Some 
of the earlier approaches have been largely discarded with the realization that blazar spectra 
are considerably more complex, characterized by an energy dependent spectral index and 
variable intensity (Samuelson et al. 1998; Aharonian et al. 1999, 2002a; Djannati-Ata'i et 
al. 1999; Krennrich et al. 1999, 2001, 2002). Construction of the intrinsic blazar spectrum 
requires therefore a detailed knowledge of the 7-ray production mechanism. For example, 
Coppi & Aharonian (1999) suggested using contemporaneous X-ray /7-ray variations in the 
blazar emission in conjuncture with a synchrotron self-Compton emission model to estimate 
the intrinsic TeV spectrum of the blazar from its X-ray spectrum. Differences between this 
intrinsic spectrum and the observed can then be used to estimate the amount of intergalactic 
absorption and the intensity of the EBL. However, there still is no consensus on the 7-ray 
emission model, which prompted an alternate approach of using either theoretical EBL 
models, or observed EBL limits and detections in order to place constraints on the intrinsic 
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blazar spectra (Guy et al. 2000; Renault et al. 2001; Dwek & de Jager 2001; Primack, 
Somerville, Bullock, & Devriendt 2001; Kneiske, Mannheim, & Hartmann 2002; Kneiske, 
Bretz, Mannheim, & Hartmann 2004; and Wright 2004 ). 

In this paper we use the most recent data on the EBL to construct twelve different 
realizations, representing all possible permutations between EBL limits and the detections 
in the different wavelength regions, to explore the effects of the EBL on the inferred spectra 
of blazars. We focus on the determination of the intrinsic energy spectra of the two most 
prominent and nearby TeV blazars, Mrk 421 and Mrk 501 which have been measured with 
good statistical precision at energies between 260 GeV and 20 TeV (Samuelson et al. 1998; 
Aharonian et al. 1999, 2002a; Krennrich et al. 1999, 2001, 2002; Djannati-Atai et al. 
1999). Some claimed EBL detections, notably the 60 /xm intensity derived by Finkbeiner, 
Davis, & Schlegel (2000), produce a pile- up of photons in the highest energy bins of the 
intrinsic blazar spectrum. This pile-up has led to claims of the existence of a "TeV-IR 
Background crisis" (Protheroe & Meyer 2000). However, the EBL spectrum required for 
calculating the 7-ray opacity towards these sources, consist largely of upper and lower limits 
and few 3<r level detections in the 0.1 - 200 /mi wavelength region (Hauser & Dwek 2001). 
Uncertainties in the EBL spectrum propagate exponentially into the derived blazar spectrum, 
a fact that has been largely ignored by previous investigators (e.g. Konopelko et al. 2003). 
Consequently, any statements claiming the existence of a "TeV-IR Background crisis" are 
unfounded considering the uncertainties in the EBL intensity and the observed 7-ray fluxes 
in the highest energy bins. 

In spite of the uncertainties in the detailed spectral behavior of the EBL, its general 
double-peak nature is well established. The peaks at ~ 1 and ~ 100 /im, are associated with 
energy releases from stars and dust, respectively. The drop in the EBL spectrum between 
these two peaks (A ~ 5 — 60 /zm) moderates the rise in the TeV opacity with 7— ray energies, 
creating an opportunity for extracting the intrinsic 200 GeV to 10 TeV spectrum of these 
nearby blazars. 

Using the observed detections and limits on the EBL we create a family of EBL spectra 
representing different realizations of the EBL, bounding extreme limits in the UV optical 
and IR-submillimeter wavelength regions. The EBL spectra used in our analysis cover a 
significantly wider range of intensities than the two "extreme" spectra considered by de Jager, 
and Stecker (2002), Konopelko et al. (2003), and the range of EBL spectra recently derived 
by Kneiske et al. (2004). The family of EBL spectral templates are used to derive the 7-ray 
optical depth to the observed blazars, to construct the intrinsic spectra of Mrk 421, Mrk 501, 
and H1426+428, and to examine spectral trends in the intrinsic variability of Mrk 421. We 
show that several intrinsic spectra display "unphysical" characteristics, manifested as an 
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exponential rise with energy, that can be directly related to the spectrum of the EBL used for 
calculating the 7-ray attenuation. Aharonian, Timokhin, & Plyasheshnikov (2002b, ATP02) 
suggested that a pile-up of photons in the highest energy bins actually represents the intrinsic 
blazar spectrum, and presented a model in which the ~ 10-20 TeV photons are produced by 
the comptonization of ambient optical photons by an ultrarelativistic jet. However, the 
exponential rises produced by some EBL realizations are unbounded, extending well beyond 
the highest energy bins in the blazar spectrum, and are therefore unlikely to be generated 
by the proposed ultrarelativistic jet-model. 

The exponential rises in the blazar spectrum derived here result from correcting the 
observed 7-ray fluxes in the highest energy bins for EBL absorption. Uncertainties in the 
observations may affect the "reality" of some of these exponential rises. We therefore employ 
the statistical F-test to examine the significance of these rises. Blazar spectra that are 
confirmed to be unphysical are then used to set new limits on the EBL. We derive the peak 
energies in the well-behaved intrinsic blazar spectra, and discuss the relevance of these results 
for blazar unification models (Costamante et al. 2001). 

Another source of uncertainty in the blazar spectrum is the absolute energy of the 7- 
ray photons. We therefore examine the effect of these uncertainties on the intrinsic blazar 
spectra and on the ability to constrain the EBL. Finally, we also examine whether the TeV 
optical depth towards Mrk 421 and 501, which are located at an almost identical redshift of 
0.03, may differ due to statistical fluctuations in the number of sources along their line of 
sight. 

The paper is organized as follows: In §2 we present the observed 7— ray spectrum of 
Mrk 421, Mrk 501, and H1426+428 and the spectral variability of Mrk 421. The current 
status of the detection and limits of the EBL are described in §3, where we also present the 
different spectral realizations of the EBL. In §4 we derive the 7-ray opacity for all the EBL 
realizations. The attenuation corrected spectra of Mrk 421 and Mrk 501 are presented in §5, 
and EBL realizations leading to unphysical intrinsic spectra are identified. The constraints 
on the EBL derived from analysis of the intrinsic blazar spectrum are discussed in §6. Blazar 
spectra corrected for absorption by the viable EBL realizations are presented in §7. A brief 
summary of the paper is presented in §8. 



2. GAMMA-RAY DATA OF MRK 421, MRK 501 AND H1426+428 

The 7-ray data used for this study were obtained by two experimental groups, the VER- 
ITAS collaboration using a single 10 m diameter imaging atmospheric Cerenkov telescope 
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at the Whipple observatory (Finley et al. 2001) in southern Arizona and the HEGRA col- 
laboration using an array of smaller size imaging atmospheric Cerenkov telescopes (Daum 
et al. 1997) on La Palma, Canary Islands. The atmospheric Cerenkov technique is prone to 
uncertainties in the absolute energy scale of the 7-ray photons resulting from systematic un- 
certainties in modeling the Earth's atmosphere, and uncertainties in the temporal variations 
of its transparency. These uncertainties limit the accuracy of the absolute energy calibration 
of the ~ TeV 7-rays to 15%, and can have an important effect on the correction applied 
to the observed blazar spectrum due to absorption by the EBL. A consistency check of the 
magnitude of major systematic uncertainties of different atmospheric Cerenkov experiments 
is given by measurements of the Crab Nebula, the standard candle of 7-ray astronomy. The 
energy spectra of the Crab Nebula measured by the Whipple telescope and the HEGRA 
detector are in agreement within systematic errors (for further details see also Aharonian et 
al. 2000; Hillas et al. 1998). 

This study is primarily concerned with the two most prominent blazars in TeV gamma- 
ray astronomy, Mrk 421 and Mrk 501. Both received a lot of attention because of their 
episodes of strong flaring activity. Furthermore, at redshifts of 0.031 and 0.034, respectively, 
they are approximately at the same distance from the observer, suggesting that absorption 
by the EBL should affect them similarly. Mrk 501 was discovered as a 7-ray blazar at TeV 
energies (Quinn et al. 1996) and to a great surprise underwent a strong and long-lasting 
flaring state during the first half of 1997 (Catanese et al. 1997; Aharonian et al. 1997; 
Djannati-Ata'i et al. 1999; Protheroe et al. 1998). The good statistics of these data led to 
the discovery of a first spectral feature in a TeV blazar, the detection of a cutoff manifested as 
a strong deviation of dN/dE from an extrapolated E~ 2 power law, at approximately 4 TeV 
(Samuelson et al. 1998). Detailed spectral measurements by other groups independently 
confirmed this cutoff feature (Aharonian et al. 1997, 1999; Djannati-Ata'i et al. 1999). The 
combined Whipple and HEGRA data cover an energy range between 260 GeV and 20 TeV, 
making the energy spectrum of Mrk 501 an interesting constraint to the EBL density. 

The second blazar that we discuss in detail in this paper is Mrk 421, which was es- 
tablished as a TeV source by Punch et al. (1992). Although strong and short flares were 
observed in 1995 and in 1996 (Gaidos et al. 1996), those data were not sufficient to unravel 
spectral features beyond a power law index (Zweerink et al. 1997; Krennrich et al. 1999). 
This changed when an active episode of long-lasting strong flares occurred for Mrk 421 in 
2000/2001, providing the best statistics ever recorded for a 7-ray blazar and also showing a 
cutoff at approximately 4 TeV (Krennrich et al. 2001; Aharonian et al. 2002a). The flaring 
data from the observing season of 2000/2001 have also been used to study the flux spectral 
index as a function of the flux level. The data have been binned according to flux levels 
resulting into 6 data sets at different flux levels (Krennrich et al. 2002). These data have 
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unraveled for the first time significant spectral variability at TeV energies. 

The spectral data of Mrk 421 and Mrk 501 from the two experiments are complementary. 
The Whipple 10 m telescope data start at 260 GeV extending out to 17 TeV for Mrk 421 
and out to 12.1 TeV for Mrk 501, whereas the HEGRA observations start at 560 GeV 
and extend to 17 TeV for Mrk 421 and to 21 TeV for Mrk 501. Due of the wealth of 
statistics, the uncertainties in the data are dominated by systematic uncertainties in the 
absolute calibration of the Whipple and HEGRA instruments. Consequently, we refrain 
from combining the data and analyze them separately. 

Finally, we present spectral data from a first more distant (,2=0.129) TeV blazar, 
H1426+428, which was discovered as a TeV source by Horan et al. (2000, 2001a, 2001b, 
2002) and confirmed by Aharonian et al. (2002a) and Djannati-Atai et al. (2002). We use 
the spectral data from the HEGRA collaboration (Aharonian et al. 2002a) and the Whipple 
collaboration (Petry et al. 2002). Because of the limited statistics of the spectra from both 
groups (5.8 and 5.5 sigma detections) we combine the spectra (as was presented in Petry et 
al. 2002). This is possible since the statistical errors dominate over the systematic uncer- 
tainties for this data set. Figure 1 summarizes the observed spectra of the three blazars, and 
presents the observed spectral variability of Mrk 421. 



3. REALIZATIONS OF THE EXTRAGALACTIC 
BACKGROUND LIGHT 

3.1. Observational Limits and Detections of the EBL 

The EBL is viewed through strong sources of foreground emission consisting of reflected 
sunlight and thermal emission from zodiacal dust particles, resolved and diffuse Galactic 
starlight, diffuse emission from interstellar dust, and emission from resolved Galactic H II 
regions and nearby galaxies (LMC, SMC). The relative contribution of these different fore- 
ground emission components varies significantly over the 0.1 to 1000 pm wavelength range 
of the EBL. 

Removal of the zodiacal dust emission from the observed sky maps posed the greatest 
challenge for the detection of the EBL in the 1 to 1000 fim wavelength region. Thermal 
emission from the zodiacal cloud dominates the foreground in the ~ 5-100 fim wavelength 
region with an intensity of vl v ~ 4 x 10 3 nW m~ 2 sr _1 between ~ 15 and 30 /im in the 
direction of the Lockman Hole. This intensity is about a factor of ~ 100 — 10 3 larger 
than the average EBL intensity expected from nucleosynthesis arguments, illustrating the 
difficulty of separating this component from the EBL. The subtraction of the zodiacal dust 



emission was described in detail by Kelsall et al. (1998). The procedure modeled the 
variation in the sky intensity caused by the Earth's motion through the IPD cloud and the 
DIRBE scanning pattern. Any zodiacal dust model is therefore insensitive to any isotropic 
emission component of the cloud. The uncertainties in the intensity of this component were 
determined from the variance in its value, obtained by modeling the primary dust cloud with 
different geometrical configurations which produced about equally good fits to the observed 
variations in the sky brightness. The Kelsall et al. (1998) zodiacal light model (KZL model 
hereafter) succeeded in subtracting 98% of the thermal emission from the zodi cloud. The 
residual 5 to 100 /im emission exhibited a strong peak at 15-30 /im. It had large systematic 
uncertainties and was far from isotropic, preventing its identification as a component of the 
EBL. 

Arguing that the ~ 15 to 30 /im residual must be mostly of solar system origin (otherwise 
the observed TeV blazars would not be detectable), Gorjian et al. (2000) used a zodiacal 
dust model similar to the KZL model, but imposed a "very strong no- zodi principle" (Wright 
2001), requiring that the 25 /im residual after the removal of the zodiacal emission be zero 
at high Galactic latitudes. We will hereafter refer to this zodiacal light model as the GWZ 
model. The resulting contribution of the GWZ model to the foreground emission is therefore 
larger than that of the KZL model, resulting in lower values for the EBL intensity in the 
1.25 to 5 /im wavelength region. 

Galactic starlight is an important contributor to the foreground emission at near infrared 
wavelengths (A ~ 1-5 /im), and the removal of this component from the DIRBE skymaps 
was discussed in detail by Arendt et al. (1998). The systematic uncertainties in the ~ 1 - 
5 /im residuals were dominated by uncertainties in the model used to subtract the emission 
from unresolved stars. Since then significant efforts have been undertaken to improve the 
removal of the Galactic stellar emission component, resulting in the detection of the EBL at 
1.25, 2.2, and 3.5 /im (Dwek & Arendt 1998, Wright & Reese 2000, Wright 2001, Cambresy 
et al. 2001, Arendt & Dwek 2003). Larger values for the EBL at these wavelengths were 
obtained by Cambresy et al. (2001) and Matsumoto et al. (2000), who used the KZL-model 
to subtract the zodiacal foreground. Use of the GWZ-model, which is characterized by a 
larger contribution of the zodiacal dust cloud to the foreground emission, will give rise to 
a lower EBL intensity. The dependence of the near-IR EBL intensity on the model used 
for the subtraction of the zodiacal emission is discussed by Arendt & Dwek (2003) and 
illustrated in Figure 6 of their paper. All COBE/DIRBE derived EBL intensities at between 
~ 1 and 2 /im are significantly higher than those derived by Madau & Pozzetti (2000) from 
galaxy number counts. This should not be too surprising, since the integrated light from 
galaxies provides only a strict lower limit on the EBL intensity, even when the integrated 
light from the resolved galaxies seems to have converged. Bernstein, Freedman, & Madore 
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(2002) pointed out that a significant fraction of the flux from resolved galaxies can remain 
undetected, since the overlapping wings of these galaxies can form a truly diffuse background 
which will be missed in standard galaxy aperture photometry. 

At mid-IR wavelengths the EBL has a lower limit at 15 /im, derived from galaxy number 
counts obtained with the Infrared Space Observatory (ISO) satellite (Elbaz et al. 2002; Met- 
calfe et al. 2003). The integrated flux from 24 /im sources detected in deep surveys with the 
Spitzer satellite gives a lower limit of 1.9±0.6 nW m~ 2 on the EBL intensity, which extrapo- 
lated to fainter flux densities provides an estimated EBL intensity of 2.7IJ7 n ^ m ~ 2 sr_1 a ^ 
that wavelength (Papovich et al. 2004). An upper limit of about 5 nW m~ 2 sr _1 was derived 
by Stanev & Franceschini (1998), Renault et al. (2001), and Dwek & de Jager (2001) from 
TeV 7-ray considerations. At far-infrared wavelengths, the EBL has been detected at ~ 200 
to 1000 /im by Puget et al. (1996), Fixsen et al. (1998), at 140 and 240 /im by Hauser et al. 
(1998), and at 100 /im by Lagache et al. (2000). Hauser et al. (1998) reported the 140 /im 
EBL intensity derived using the DIRBE photometric calibration. A somewhat lower value 
(but consistent with the DIRBE calibration) is derived if the FIRAS photometric scale is 
used in the calibration (Hauser et al. 1998, Hauser & Dwek 2001). 

At UV wavelengths Gardner et al. (2000) presented the integrated light obtained from 
deep galaxy counts using the Space Telescope Imaging Spectrograph (STIS) combined with 
counts obtained with the FOCA balloon-borne telescope (Milliard et al. 1992). Bernstein, 
Freedman, & Madore (2002) reported the detection of the EBL at 0.3, 0.55, and 0.8 /im, at 
levels that are higher than the results of Madau & Pozzetti (2000), results that were disputed 
by Matilla (2003). Current limits and detections of the EBL are presented in Table 1, and 
depicted in Figure 2. 



3.2. Template EBL Spectrum 

Given the general double-peak nature of the EBL and the uncertainties in its spectral 
intensity we created a family of EBL spectra spanning extreme possible combinations of rel- 
ative peak values. Altogether we constructed twelve template spectra, representing different 
realizations of the EBL, by fitting polynomials to all possible combinations of the following 
spectral components: 

1. Three stellar components consisting of: (1) high-UV: defined by the 0.1595 and 
0.2365 /im data from Gardner et al. (2000), and the Cambresy et al. (2001) de- 
terminations that used the KZL zodi-light model; (2) mid-UV: defined by the 0.1595 
and 0.2365 /tm data from Gardner et al. (2000), and the Wright (2001) EBL intensity 
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derived by using the GWZ zodi-light model; and (3) low-UV: denned by the 0.1595 
and 0.2365 /im data from Gardner et al. (2000), and the number counts of Madau 
& Pozzetti (2000). The three realizations of the stellar components are driven by 
the near-IR and the 0.1595 and 0.2365 /im constraints, and are compatible with the 
claimed detections by Bernstein et al. (2002), considering the uncertainties in their 
values. 

2. Two mid-IR components, which were defined by the uncertainties in the 15 /im lower 
limit determined by Elbaz et al. (2002) from ISO observations. The high mid-IR 
EBL is represented by the nominal 15 /im intensity + 3<r, and the low mid-IR EBL is 
represented by the nominal 15 /im intensity — 3a. 

3. Two far-IR components, defined by the 2 different calibration of the DIRBE 100 and 
240 /im data points. The DIRBE calibration gives rise to higher values of the EBL 
compared to the FIRAS calibration (see Table 1). Above ~ 240 /im all components 
were fitted to the FIRAS determination of Fixsen et al. (1998). 

Finally, we also constructed an "average" EBL spectrum, defined by a polynomial fit through 
the nominal 15 /iin lower limit, and the average UV and far-IR limits. It is approximate in 
nature, and only used for illustrative purposes. 

For sake of brevity the twelve different EBL realizations will be referred to as XYZ, with 
X=L, M, or H representing the intensity (low, medium or high) of the stellar component of 
the EBL, Y=L or H, representing the low or high intensity of the 15 /im EBL flux, and Z=L 
or H, representing the low or high intensity of the EBL flux at far-IR wavelengths. So an 
EBL realization designated as MHL represents an EBL derived by a polynomial fit through 
the mid-UV, the high mid-IR, and the low far-IR spectral intensities of the EBL. 

Figure 2 shows the different realizations of the EBL, and the observational data points 
used in their derivation. The three components are represented in the figure by solid lines 
(high-UV), lines with connected dots (mid-UV), and dashed lines (low-UV). Lines going 
through the high mid-IR are black, whereas those going through the low mid-IR point are 
grey. All lines going through the high far-IR data are thick compared to the ones going 
through the low far-IR data points. Also shown in the figure is a nominal "average" EBL 
spectrum which is represented by a heavy dark line going through the 15 /iin lower limit. For 
sake of comparison we also show the range of EBL intensities sampled by the two EBL spectra 
chosen by de Jager & Stecker (2002) and Konopelko et al. (2003) to represent the EBL limits 
(shaded area in top figure). The mid- and far-IR intensities of these spectra were derived 
from simple backward evolution models using highly idealized galaxy spectra (Malkan & 
Stecker 2001). At UV and optical wavelengths, de Jager & Stecker (2002) augmented these 
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spectra with the integrated galactic light derived by Madau & Pozetti (2000). As is evident 
from the figure, these two spectra sample a very small range of viable EBL spectra. The 
bottom panel compares our different EBL realizations with the range of EBL spectra derived 
by Kneiske et al. (2004) for various cosmic star formation histories. Their models do not 
reproduce the high EBL intensities at UV to mid-IR wavelengths, a common problem in all 
EBL models. The data points depicted in the figure are listed in Table 1, and the coefficients 
of the polynomial approximations are given in Table 2. 



4. THE 7-RAY OPACITY OF THE LOCAL UNIVERSE 

The cross section for the 7 + 7 — > e + + e~ interaction of a 7-ray photon of energy E 1 
emitted from a source at redshift z with a background photon of energy e is given by (e.g. 
Jauch & Rohrlich 1955) 
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where <jt = 6.65 x 10~ 25 cm 2 is the Thompson cross section, e t h the threshold energy of the 
interaction, and \x = cos 9, where 9 is the angle between the incident photons. The 7-7 cross 
section for the interaction with an isotropic distribution of background photons has a peak 
value of 1.70xl0~ 25 cm 2 for /3 = 0.70, which corresponds to energies for which the product 
E^e p» 4(m e c 2 ) 2 « 1 MeV 2 , or A e (/im) m 1.24£' 7 (TeV), where X e is the wavelength of the 
background photon. 

The optical depth traversed by a photon observed at energy E 1 that was emitted by a 
source at redshift z is given by: 
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where n e (e f , z')de' is the comoving number density of EBL photons with energies between 
e' and e'+de' at redshift z', e' th = tthiE'^-, p), E' = E 7 (l + z'), and where d£/dz, is given by 
(e.g. Peacock 1999): 
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e(z) = {(i + z) 2 (n m z + i) + z(2 + z)[(i + z ) 2 n r -n A ]} 1/2 

where Q m and Q r are, respectively, the matter and radiation energy density normalized to 
the critical density, Q A = A/3H^ is the dimensionless cosmological constant (fi m + Q r + Qa 
= 1 in a flat universe), Rh = c/H is the Hubble radius, c is the speed of light, and H is 
the Hubble constant, taken here to be 70 km s _1 Mpc -1 . The comoving number density of 
EBL photons of energy e at redshift z is given by: 

e 2 n e (e, z) = I — ) vl v {y, z) (4) 



i/CJi/, z'\ 



dt 



dz> 



dz' 

l + z> 



where e = his, v' = v{\ + z'), and C v {v' ', z') is the specific comoving luminosity density at 
frequency i/ and redshift z 1 . 

Figure 3 (right panel) depicts the TeV opacity of a source located at redshift z = 0.030 
to background photons with an EBL spectrum given by the average spectrum depicted in 
the left panel of the figure (see also Figure 2). The shaded curves in the figure represent the 
contribution of the different wavelength regions (depicted in the shaded bar diagram in the 
left panel) to the total opacity. The figure illustrates the relation between the rate of increase 
in the TeV opacity with 7-ray energy, and the spectral behavior of the EBL. Particularly 
noticable is the decrease in the rise of the opacity between ~ 1 and 5 TeV, resulting from 
the dip in the EBL intensity between the stellar and dust emission peaks. 

As an aside, we note that in calculating the 7-ray opacity we discovered a numerical 
error of about 40% in the polynomial approximations presented by de Jager & Stecker (2002) 
which apparently was caused by the coarse grid used in the integration of r 77 over angles 
(de Jager, private communication). 



5. NEW GAMMA-RAY DERIVED LIMITS ON THE 
EXTRAGALACTIC BACKGROUND LIGHT INTENSITY 

The spectral energy distribution of blazars consists of two spectral components: (1) a 
low-energy X-ray component, extending up to energies of about 100 keV which is attributed 
to synchrotron radiation from energetic electrons, and (2) a high-energy 7-ray component 
with energies extending to the TeV range, which is usually attributed to inverse Compton 
(IC) scattering of the synchrotron emission by energetic electrons (Maraschi, Ghisellini & 
Celotti 1992; Marscher & Travis 1996). Competing models exist for the nature of the particles 
producing the 7-ray emission. Leptonic models assume that energetic electrons are the 



-13- 



primary particles producing the 7-rays by IC scattering off the synchrotron emission or 
other ambient soft photons. Hadronic models assume that the 7-ray emission is produced 
by proton-induced synchrotron cascades or by decays of secondary particles such as neutral 
pions and neutrons (Mannheim & Biermann 1992; Mannheim 1993, 1998) or, alternatively, 
by synchrotron radiating protons (Miicke & Protheroe 2001; Aharonian 2000). 

All the above models predict a decline in the blazar luminosity at the highest 7-ray 
energies. Specifically, none of these models predict an exponential rise of luminosity with 
energy. An exception is the model of Aharonian, Timokhin, & Plyasheshnikov (2002b) which 
was constructed to explain a possible pile-up in the highest energy bins of the intrinsic blazar 
spectrum. In their model the intrinsic blazar spectrum rises sharply between ~10 and 20 
TeV and has an abrupt cutoff at ~20 TeV. Coppi (2003) emphasizes that such an upturn in 
luminosity at TeV energies is not impossible but unlikely, since it would require extremely 
energetic particles in the jet. Furthermore, the exponential rises produced by some EBL 
realizations are unbounded and therefore unlikely to be generated by the jet model. Since 
most blazar models are moderately successful in explaining the gross features of the low and 
high energy emission peaks in the blazar spectra over a wide range of energies, we consider 
any exponential rise in a blazar spectrum at TeV energies as "unphysical" . 



The intrinsic photon spectrum, (dN/dE)i, of a blazar located at redshift z is given by: 

g^W^lx^ (5) 

where (dN/dE) b s , the observed spectrum, can be fit with a function f b s consisting of a 
power law with an exponential cutoff of the form: 

f obs = <5>E~ a x exp[-E/E ] (6) 

Using the spectral templates representing the various realizations of the EBL described 
in §3, we calculated the 7-7 opacity towards Mrk 421 and 501, and derived their intrinsic 
7-ray spectra. The intrinsic blazar spectrum can be fit by either a parabolic function, f p , or 
a parabolic function with an exponential rise, f e , of the form: 



dE 



f p (E, $, a, (3,) = ^E- a ~^^E) (7) 

f e (E, $, a, [3, E Q ) = $ E- a -^ lo ^(E) x exp(£/£ ) (8) 



The function f e is a generalization of f p , designed to explore intrinsic spectra that show an 
"unphysical" behavior, characterized by a sudden exponential rise in the function E 2 (dN/dE)i 
after an initial monotonic decline or flat behavior with energy. 
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We examined the statistical significance of an exponential rise in the intrinsic blazar 
spectrum by using the F-test to calculate the probability that the reduction in the x 2 of 
the fit due to the inclusion of the additional parameter E exceeds the value which can 
be attributed to random fluctuations in the data. If the reduction in the x 2 is sufficiently 
large then the exponential rise is statistically significant, the intrinsic spectrum is considered 
unphysical, and the EBL spectrum causing this behavior is excluded as a viable realization 
of the EBL. 

Let F x be the ratio Ax 2 /xl 2 = [xVi) - X \^)} / 'xl 2 , where xl x and xl 2 are the re- 
duced x 2 for i>i = N — m and v 2 = N — (m + 1) degrees of freedom, respectively. The 
value of F x measures the fractional improvement in x 2 , and its statistics follow that of the 



F(i/i, v 2 ) = (Xi/^iVCx!/^) distribution with v\ = 1, and u 2 = N — (m + 1). Each F, 



x 



value has therefore an integral probability distribution P(F X , v\ = 1, v 2 \ which measures 
the probability that the improvement in the fit was not a random event. Table 3 list the 
resulting probabilities for the F-test on F x . A P-value of 95% or larger is commonly regarded 
as significant. Consequently, EBL realizations giving rise to intrinsic spectra for which the 
exponential rise is significant, are ruled out as viable spectral representations of the EBL, and 
are labeled "0" in the table. Acceptable EBLs are designated with a "1". The table shows 
that Mrk 501 provided more stringent constraints on EBL scenarios. This is not surprising 
considering the fact that the 7-ray observations of this blazar extend to higher energies and 
have smaller error bars, compared to the observations of Mrk 421. 

Figures 4 and 5 depict the intrinsic spectra of Mrk 421 and 501, respectively, for all 
EBL realizations. The figures clearly show the "unphysical" behavior of the blazar spectra 
for the LHH, MHH, MHL, and HHL realizations of the EBL. For both sources the intrinsic 
spectra for these realizations initially decrease with energy, then suddenly increase for at 
least two data points at energies above ~ 10 TeV. The rise can be directly attributed to the 
EBL intensity in the wavelength regime contributing to the 7-ray opacity. All, except two, 
of the EBL realization with a high mid-IR intensity give rise to unphysical blazar spectra, 
and are therefore ruled out by the 7-ray observation. Of the two viable EBL realizations 
with high mid-IR intensity, LHL has a low UV and low far-IR intensity, and HHH has high 
UV and high far-IR intensities. So the common characteristic of these two spectra is that 
they do not connect between opposite extremes in the UV and far-IR intensities of the EBL, 
that is, high UV to low far-IR and vice versa. 

Using only visual inspection of the intrinsic blazar spectra (e.g. Figure 5) , one would 
be tempted to regard additional realizations such as HHH, LLH, MLH, and LHL unphysical, 
since they also give rise to an increasing 7-ray spectrum. However, the rises occur only at 
the last datum point, and are statistically insignificant. This fact is confirmed by the results 
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of the statistical F-test that show that the confidence level of these rising spectra is not high 
enough to ascertain their reality. These examples illustrate the importance of conducting 
rigorous statistical tests over casual visual inspections, for determining the reliability of 
features in the derived intrinsic spectrum of blazars. 

In addition to the systematic uncertainties of the individual spectral data points we also 
considered the possibility that the reported gamma-ray energies are higher than their actual 
values, within the claimed experimental uncertainties. Use of the published higher energy 
values would result in an overestimate of the gamma-ray opacity and could cause an artificial 
rise in the absorption corrected spectra. We have adopted the conservative assumption that 
the reported energy spectra of Mrk 501 and Mrk 421 are shifted up in energy by 15%. We 
therefore recreated the observed HEGRA spectra for these sources by adopting an energy 
grid that is 0.85 times the reported one, and derived our constraints for the EBL for these 
spectra using the previously described F-test. The results show that three previously rejected 
EBL realizations: LHH, MHL, and HHL are now viable, and only the MHH one can be firmly 
rejected. 

Figure 6 depicts the effect of the energy shift on the intrinsic spectrum of Mrk 501. The 
open squares represent the nominal blazar spectrum, and the solid ones represent the blazar 
spectrum when the photon energies are shifted down by 15%. The heavy dashed and solid 
lines represent functional fits (a power-law with an exponential cutoff, eq. 6) through the 
data. The intrinsic spectra for both cases are calculated using the LHH realization of the 
EBL and are, depicted by open and solid circles for the nominal and energy-shifted blazar 
spectrum, respectively. The dashed line is a functional fit through the intrinsic spectra of 
the blazar using its nominal observed spectrum (see also Fig. 5a). As mentioned before, 
the exponential rise is statistically significant, and the LHH realization of the EBL was 
rejected. The two solid lines represent parabolic and exponentially rising fits to the intrinsic 
blazar spectrum, derived from its energy shifted spectrum. The figure visually shows that 
the exponentially rising function does not provide a statistically better fit to the intrinsic 
blazar spectrum than the parabolic one. Hence, the LHH realization of the EBL cannot 
be rejected in this case. These results clearly demonstrate that the absolute calibration of 
atmospheric Cherenkov telescopes directly impacts the ability to constrain the EBL with 
7-ray astronomy. It is therefore important that systematic uncertainties be improved in the 
next generation of telescopes. 

Figure 7 depicts the eight viable EBL spectra (HHH, LLH, MLH, HLH, LHL, LLL, MLL, 
and HLL) that give rise to physically "well behaved" intrinsic 7-ray spectra. In addition 
we show as dashed lines, the three EBL realizations (LHH, MHL, and HHL) that are viable 
when the observed 7-ray energies are shifted down by 15%. Figure 8 shows the optical depths 
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derived for these EBL realizations. All optical depths are characterized by a flattening in 
the 1 to 10 TeV optical depth caused by the dip in the EBL at mid-IR wavelengths. This 
is in contrast to the model of Konopelko et al. (2003) which shows a monotonic increase 
of 7-ray opacity with energy for their adopted EBL spectra. The shaded area in the figure 
depicts the optical depth bounded by the two EBL models adopted in their analysis, and 
depicted in Figure 2 (upper panel) of this paper. The heavy shaded area represents the range 
of opacities corresponding to the EBL spectra of Kneiske et al. (2004) shown in Figure 2 
(lower panel). 

From Figure 8 we see that all EBL realizations imply substantial absorption at energies 
above 400 GeV for blazars with a redshift of more than z ~ 0.03 except for the LLL case 
(lowest curve in the figure) for which the spectral shape is only marginally modified by the 
attenuation. The energy dependence of the opacity is crucial in determining the shape of the 
intrinsic blazar spectrum and the position of its energy peak. This is especially important 
for modeling the spectral variability of Mrk 421, for which observations suggest that its peak 
energy is shifting to higher TeV energies with increasing flaring intensity, and for determining 
the intrinsic spectrum of H1426+428, which is located at z=0.129. 



6. IMPLICATIONS FOR THE EXTRAGALACTIC 
BACKGROUND LIGHT INTENSITY 

6.1. New Limits on the EBL Spectrum 

Using minimal assumptions on the behavior of the intrinsic blazar spectrum we derived 
new limits on the spectrum of the EBL. Figure 7 shows the eight template spectra that 
did not give rise to an "unphysical" behavior in the blazar spectrum. Table 4 summarizes 
key characteristics of the eight viable EBL realizations: the total EBL intensity over the 
0.1 to 10 3 /im wavelength range, and the intensities attributed to starlight (0.1-5 /im) and 
reradiated thermal emission from dust (5-10 3 /im). The last column gives the fractional 
contribution of the dust emission to the total EBL intensity. 

Finkbeiner, Davis, & Schlegel (2000) detected an isotropic emission component at 60 
and 100 /im after the subtraction of the expected zodiacal and Galactic emission components 
from the COBE/DIKBE maps. The intensity of the residual emission was found to be 28.1 ±7 
and 24.6±8 nW m~ 2 sr" 1 at these respective wavelengths. Our analysis puts a strong upper 
limit of ~ 15 nW m -2 sr _1 on the EBL intensity at 60 /im, otherwise the EBL will give rise 
to an unphysical blazar spectrum. This upper limit is still within the la uncertainty of their 
nominally stated detection. However, any isotropic components above this upper limit must 
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be of local origin. 

The maximum EBL intensity corresponds to the HHH realization is 140 nW m -2 sr -1 , 
and the minimum EBL intensity, corresponding to LLL, is 48 nW m~ 2 sr _1 . The HLL 
realization gives a fractional dust contribution to Itot of only 23%, significantly less than the 
fractional contribution of dust to the local luminosity density, which is about 30% (Soifer & 
Neugebauer 1991, Dwek et al. 1998). The fractional contribution of the dust emission to the 
EBL intensity should be larger than its contribution to the local luminosity density, since in 
the past galaxies emitted a larger fraction of their total luminosity at infrared wavelengths, 
as suggested from the evolution in the number counts of ultra-luminous IR galaxies with 
redshift (e.g. Chary & Elbaz 2001). So the HLL realization of the EBL may be rejected 
on spectral grounds. All other seven EBL realizations give I dust I 'hot values larger than 0.34, 
which is marginally larger than the local value. 



6.2. Local Fluctuations in the Number Density of Background Photons 

The two blazars Mrk 421 and 501 have almost identical redshifts (z = 0.031 and 0.034, 
respectively), but they may have different opacities if local fluctuations in the number of 
infrared galaxies are important. The optical depth to 7-ray photons going through an indi- 
vidual disk galaxy with a mid-IR luminosity Lj R , mid-IR photon number density rim, and 
disk radius R is approximately given by: 

Tgai ~ o" 77 n IR R (9) 

fl L IR /e IR ^ 




vC R 2 
- 0.2 a 25 L 10 R- p \ 

where (T25 is the 7-7 cross section in units of 10~ 25 cm 2 , L w is the mid-IR luminosity of the 
galaxy in units of 10 10 L Q , Rk pc is a typical galactic radius (in kpc) containing most of the 
IR emission from the galaxy, and em is the energy of a mid-IR photon with a wavelength of 
10 /im. Equation (9) shows that for an appropriate combination of Lj R and R, individual 
galaxies can be optically thick to TeV 7-rays. For example, the optical depth will be unity 
for galaxies with a mid-IR luminosity of ~10 12 L and a disk radius of ~4 kpc. So a single 
galaxy along the line of sight to a blazar can have a significant affect on its observed spectrum. 
However, the probability of encountering such galaxy along the line of sight to either Mrk 421 
or Mrk 501 is very small. The local IR luminosity density is about 2 x 10 7 L Q Mpc~ 3 (Dwek 
et al. 1998). Assuming that all this energy is emitted at mid-IR wavelengths gives an 
upper limit of n ga i ~ 10~ 5 Mpc~ 3 on the number density of these TeV-thick galaxies. The 
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probability for intersecting such galaxies is given by P fa n ga iR 2 L, where L fa 130 Mpc is 
the distance to z — 0.03. The resulting probability for encountering a galaxy opaque to TeV 
7-rays is about 2xl0 -8 . It is therefore statistically very unlikely that random fluctuations 
in the number of galaxies along the line of sight to these two TeV blazars will cause any 
significant differences in their TeV opacity. 



7. THE ABSORPTION-CORRECTED BLAZAR SPECTRA 

7.1. The Intrinsic Average Spectrum of Mrk 421 and Mrk 501 

Figure 9 shows the absorption-corrected spectra of Mrk 421 applied to the Whipple 
and HEGRA y-ray data using the eight viable EBL realizations depicted in Figure 7. The 
dotted and dashed lines represent the intrinsic blazar spectra derived by using, respectively, 
the HLL and LHL realization of the EBL. The former was rejected on the basis of spectral 
considerations: the fraction of the total EBL intensity radiated by dust was too small and 
inconsistent with the fact that the infrared luminosity density increases with redshift. The 
latter, as can be seen in Figure 5g, gives rise to an unphysical blazar spectrum, although 
this unphysical behavior was found to be statistically insignificant. 

We note, that all absorption-corrected Mrk 421 spectra exhibit a curved shape. Assum- 
ing that there is no additional attenuation by radiation fields in the proximity of the blazar 
jet, these represent the range of possible intrinsic spectra, from which we can derive their 
peak luminosity, which provides important constraints for models for the origin of the y-ray 
emission. The figures show that the luminosity peak of Mrk 421 occurs within the energy 
range between 0.5 TeV and 1 TeV. 

Mrk 501 is approximately at the same redshift as Mrk 421, hence EBL absorption should 
affect its spectrum similarly, allowing a direct comparison of the intrinsic spectra of these 
two blazars. Figure 10 shows absorption-corrected spectra of Mrk 501 for acceptable EBL 
realizations. The dotted and dashed line represent the intrinsic spectra for the same EBL 
realizations as in Figure 9. The unphysical nature of the LHL-corrected intrinsic spectrum 
is clearly shown for the Mrk 501 HEGRA data: the intrinsic spectrum is essentially flat with 
energy. The results from HEGRA and Whipple data show that the intrinsic averaged energy 
spectrum of Mrk 501 during its flaring activity in 1997 peaked at energies between 0.8 and 
2.5 TeV. 

A more rigorous determination of the peak energy and its uncertainty can be made from 
the parabolic fit to the spectra of these two blazars (equation 7). Since the uncertainties 
in the derived peak energies are dominated by the wide range of EBL realizations, and 
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no truly contemporaneous X-ray data are available for the complete data set, this simple 
function is sufficiently accurate to get an estimate of the peak energy of the blazar. The use 
of more elaborate fitting functions based on detailed blazar models may be justified once 
data obtained with the next generation of atmospheric Cerenkov telescopes covering a wider 
range of 7-ray energies and contemporaneous X-ray observations will become available. 

The peak energies, E pea k, for Mrk 421 and 501 are presented in Table 5. For Mrk 421 
both, the HEGRA and the Whipple data sets independently show that the intrinsic time- 
averaged spectrum of Mrk 421 during the 2000/2001 season peaked at an energy between 
0.5 and 1.2 TeV. For Mrk 501, the HEGRA data give a consistently higher peak energy than 
the Whipple data, ranging from 0.7 to 1.8 TeV for the Whipple, and from 1.1 to 2.5 TeV 
for the HEGRA data. Furthermore, the table shows that peak energies of Mrk 501 are 
systematically higher by about 30 to 60% compared to that of Mrk 421 regardless of the 
EBL realization used to derive their intrinsic spectra. This point is also illustrated in Figure 
11 which depicts the Mrk 501 versus the Mrk 421 peak energies for all acceptable EBL 
realizations. Table 6 shows the peak energy in the intrinsic spectrum of Mrk 501 obtained 
by correcting the energy-shifted spectrum for absorption by the different EBL realizations. 
Comparison to the peak energies in Table 5, shows that the general trend is to shift the peak 
energies to lower values. Uncertainties in the determination of the absolute photon energy 
therefore affects the intrinsic spectrum and the determination of peak energies. 

The small difference of < 60% in the peak energies of these two blazars while they were 
in a high flaring state is quite surprising considering the fact that their historical synchrotron 
spectra in X-rays have peaked at largely different energies: at E > 100 keV for Mrk 501 (Pian 
et al. 1998; Catanese et al. 1997), and at E ~ 0.5-8 keV for Mrk 421 (Brinkmann et al. 2000; 
Krawzcynski et al. 2001; Tanihata et al. 2004). The measurements of the synchrotron peak 
energy of Mrk 501 were carried out during the same months as the 7-ray data used in this 
paper were recorded (see also Catanese et al. 1997; Samuelson et al. 1998). However, the 
published peak synchrotron energy of Mrk 421 was not measured contemporaneously with the 
7-ray data. Spectral measurements at X-ray energies of Mrk 421 during its 2000/2001 flare 
do exist, but have not yet been made publicly available. Preliminary results from this flaring 
state suggest that the location of the spectral peak is consistent with those determined in 
previous years (Krawzcynski private communication), suggesting that the synchrotron peaks 
of the two blazars during their high flaring state are indeed largely different. 
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7.2. The Intrinsic Spectral Variability of Mrk 421 

The measurement of spectral variability of Mrk 421 (Krennrich et al. 2002) allows 
further studies of its time-averaged spectra over a range of flux levels between 1 and 10 Crab. 
The observed spectra of Mrk 421 show an increased flattening of the spectra with increasing 
excitation level, suggesting a systematic shift of the peak energy to higher energies as a 
function of excitation level. To determine the peak energy of their intrinsic spectrum we 
corrected the observed spectra for the different epochs for absorption by all the eight viable 
EBL realizations derived in §5. 

Figure 12 (top left panel) shows the average energy spectra of Mrk 421 of various 
flaring states with the data binned by flux level based on measurements with the Whipple 
observatory 10 m 7-ray telescope. The remaining panels show the absorption-corrected flare 
spectra using the acceptable EBL realizations. It can be clearly seen that the EBL corrected 
spectra imply higher peak energies for the high emission states than the un-corrected data. In 
fact the peak energy for the highest flux level ranges from 730 GeV (LLL) to approximately 
2 TeV (HHL). 

Figure 13 depicts the evolution of the peak energy with flaring state. All energies are 
normalized to the peak energy of the lowest state, except for the LHL case for which the peak 
energies in the first and third flaring state could not be determined. The figure shows that 
for all EBL realizations the peak energy shifts towards higher energies with increasing flaring 
state. Table 7 lists the peak energies for the different flaring states and EBL realizations. 
The intrinsic spectrum of Mrk 421 in the lowest flaring state has a peak energy that falls 
significantly below the energy range covered by the data and therefore cannot be accurately 
determined, as evident by the large error bars. However, in spite of the uncertainty in 
the peak energy, the figure does show that for several EBL realizations, the drop in the 
peak energy of this flaring state is significant. Also noticeable in the figure is the relative 
flattening in the value of the peak energy between the second and fourth flaring state. A 
similar behavior has been observed for Mrk 421 at X-ray energies (George et al. 1988; Fossati 
et al. 2000; Brinkmann et al. 2000) showing low and soft X-ray states and hard and bright 
high emission states. 



7.3. The Intrinsic Spectrum of H1426+428 

A search for other extreme TeV blazars similar to Mrk 501 with synchrotron peak 
energies in the hard X-ray band was carried out by Costamante et al. (2001) using the 
BeppoSAX satellite. They identified the BL Lac object H1426+428 as a strong hard X-ray 
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source. Spectral measurements of H1426+428 in February 1999 using the LECS, MECS and 
PDS instruments on BeppoSAX revealed a peak energy above 100 keV, similar to Mrk 501 
during a long lasting outburst in 1997 (Catanese et al. 1997; Pian et al. 1998). This was 
followed by the first weak detection of a 7-ray signal above 300 GeV by Horan et al. (2000, 
2002). Consequently, H1426+428 has become an interesting comparison object to the TeV 
blazars Mrk 421 and Mrk 501. The third extreme blazar with a high synchrotron peak energy 
(Giommi, Padovani, & Perlman 2000) and shown to emit TeV 7-ray radiation (Catanese et 
al. 1998) is 1ES2344+514, but a spectral analysis in 7-rays is not available yet. 

H1426+428 has a redshift of z=0.129, therefore its 7-ray spectrum above 400 GeV is 
more strongly attenuated by the EBL than in the case of Mrk 421 and Mrk 501. Our results 
show that for various EBL realizations, the Mrk 421 and Mrk 501 spectra above 400 GeV 
are already substantially attenuated (see Fig. 8). A source four times more distant should 
be heavily absorbed above 400 GeV. As a consequence, a direct comparison of the intrinsic 
spectrum of H1426+428 with Mrk 421 and Mrk 501 is more complicated since it depends on 
the evolution of the EBL in the ,2=0-0.13 redshift interval. 

Having limited the number of EBL realizations using the Mrk 421 and Mrk 501 data, we 
applied the remaining viable EBL realizations to H1426+428 using the spectral data from 
Petry et al. (2002) and Aharonian et al. (2003). In our calculations we assumed that the 
EBL is constant in the to 0.13 redshift interval, so that the 7-ray opacities given in Figure 8 
were simply rescaled to the distance of H1426+428. The absorption-corrected spectra are 
shown in Figure 14, and the curves in the plots represent analytical approximations to the 
observed and intrinsic source spectra. Observed as well as intrinsic spectra were fitted by a 
power law (eq. 6 with Eq — ► 00), and are shown as solid lines. However, intrinsic spectra 
corrected for the HLH, LHL, MLL, and HLL realizations of the EBL were better fitted by a 
parabolic function (eq. 7), and are depicted by a dashed line in the figure. The figure shows 
that there is a large spread in peak energies for the different EBL realizations. The LLL and 
LLH realizations yield peak energies below 400 GeV, the MLL, HLL, MLH and HLH cases 
provide peak energies of 1-5 TeV whereas the HHH realization suggests a luminosity peak 
above 10 TeV. 

This range of possible peak energies is significantly different than the one derived for the 
7-ray spectrum of H1426+428 in previous studies. For example, Aharonian et al. (2003) use 
various theoretical models for the EBL as well as an extreme phenomenological spectrum 
to derive an either flat [in E 2 (dN/dE)] or a strongly rising intrinsic energy spectrum of 
H1426+428. Using the Primack et al. (2001) EBL models and an EBL spectrum designed 
to match the UV-optical upper limits Costamante et al. (2003) derive an intrinsic spectrum 
that peaks at energies above 8-10 TeV. Our work, using the most current observational limits 
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and detections of the EBL, suggests alternative possibilities including a strongly rising, an 
extremely flat, or a parabolic-shaped intrinsic energy spectrum in the 1-5 TeV regime. Six 
out of the eight EBL realizations give rise to a soft intrinsic spectrum or a luminosity peak 
around 1-5 TeV. Only the extreme HHH and LHL cases suggests a luminosity peak above 
10 TeV. 

In fact, three of the EBL scenarios (HLH, MLL, and HLL) yield an intrinsic energy 
spectrum for H1426+428 that exhibits a peak energy of 1-5 TeV similar to the intrinsic 
7-ray spectra of Mrk 421 and Mrk 501. This concurs well with the fact that Mrk 501 also 
has a synchrotron peak energy similar to H1426+428. 

Only the HHH and LHL realizations lead to a strongly rising intrinsic spectrum with 
energy. These same realizations also yielded an unphysical rise in the intrinsic spectrum 
of Mrk 501 (see Fig. 5), but these rises were not considered to be statistically significant 
by the F-test. Consequently, the question as to whether or not the TeV 7-ray spectrum of 
H1426+428 exhibits a dramatic rise with peak energies above 10 TeV remains unsolved. A 
rising spectrum with a peak energy above 10 TeV cannot be excluded, but is certainly not 
favored over the six other EBL realizations leading to peak energies in the few TeV regime 
or below. 

It is clear from the range of absorption-corrected spectra of H1426+428, that more de- 
tailed spectral measurements in the 7-ray regime would be necessary to resolve the peak 
energy of H1426+428 and further narrow the viable range of EBL realizations. The large 
systematic error in the absorption-corrected spectrum results from the fact that H1426+428 
is heavily absorbed above 400 GeV and uncertainties in the EBL realizations enter exponen- 
tially in the derivation of its absorption-corrected spectrum. The uncertainty in the EBL 
spectrum is the dominant factor in reconstructing the peak energy of the H1426+428 intrinsic 
spectrum. This is clearly different from the case of Mrk 421 and Mrk 501 for which we were 
still able to provide an estimate of the intrinsic peak energies, despite the large uncertainties 
in the EBL. Blazars like H1426+428 with redshifts of z > 0.1 will be extremely useful for 
constraining the EBL once measurements by GLAST, HESS, MAGIC and VERITAS will 
be available with sensitivity starting in the sub-100 GeV regime, where the transition from 
unabsorbed to the absorbed part of the 7-ray spectrum is expected to occur. In this region 
it should again be possible to limit the EBL cases by testing for unphysical rises in the 7-ray 
spectrum. 
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8. DISCUSSION AND SUMMARY 

In this paper we presented the observed TeV 7-ray fluxes from the three blazars Mrk 421, 
Mrk 501, and H1426+428, and the evolving spectrum of Mrk 421 during a period of intense 
variability. The observed blazar fluxes are presented in Figure 1. The intrinsic spectrum from 
these sources is attenuated on route to earth by low energy UV to submillimeter wavelength 
photons that constitute the extragalactic background light (EBL). The main objective of this 
paper was to explore the range of possible intrinsic source spectra of these nearest blazars, 
and derive new limits on the spectrum of the EBL, using the most recent constraints on the 
EBL spectrum, and minimal assumptions regarding what constitutes an unphysical behavior 
in the blazar spectrum. 

Using observed limits and detection of the EBL we constructed a family of twelve 
realizations of the EBL spectrum, and derived the intrinsic spectra of the three blazars. 
The different spectral templates of the EBL are presented in Figure 2. The intrinsic spectra 
of Mrk 421 and 501 derived for all the EBL realizations are depicted in Figures 4 and 5, 
respectively. Some EBL realizations led to an unphysical behavior in the intrinsic blazar 
spectrum, characterized by an exponential rise after a decline or flat behavior with energy. 
In §5 we apply the F-test to examine the statistical significance of the exponential rise in an 
intrinsic blazar spectrum. Table 3 summarizes the results of the F-test, flagging out EBL 
realizations that are considered unphysical. 

Discarding EBL realizations that lead to unphysical 7-ray spectra, leaves eight viable 
EBL spectra, shown in Figure 7. Most EBL spectra with high mid-IR intensity led to un- 
physical blazar spectra. Table 4 summarizes some characteristics of these EBL realizations. 
The range of EBL intensities is between ~ 50 and 140 nW m -2 sr _1 . The fractional contribu- 
tion of dust emission to the total EBL intensity ranges from 0.23 to 0.66. The contribution 
of the dust emission to the total luminosity density in the local universe is about 0.30. Since 
in the past galaxies reradiated a larger fraction of their stellar energy output at infrared 
wavelengths, the spectrum leading to a fractional dust contribution of 0.23 to the total EBL 
intensity can be ruled out as a viable realization of the EBL. 

Figure 8 depicts the 7-ray opacity of the universe to z — 0.03 for Hq = 70 km s _1 
Mpc^ 1 . The opacity exhibits a flattening at energies between 1 and 10 TeV, resulting from 
the dip in the EBL intensity at mid-IR wavelengths. The 7-ray opacities we derive show 
markedly different behavior with energy compared to the opacities derived by Konopelko et 
al. (2003), which are monotonically increasing with energy, reflecting the EBL spectra used 
in their calculations. 

Figures 9 and 10 depict the intrinsic spectra of Mrk 421 and 501 for each of the eight 
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remaining EBL realizations. Whereas the observed spectra exhibited a monotonic decrease 
with 7— ray energy, all absorption corrected spectra exhibit a parabolic behavior with energy. 
This allowed us to identify a peak energy, E vea ^^ in the spectra of these blazars. For all EBL 
realizations we found that the peak energy of Mrk 421 was between 0.5 and 1.2 TeV, whereas 
that of Mrk 501 was between 0.8 and 2.5 TeV, depending on the EBL realization. The value 
of E pea k of Mrk 501 was consistently higher than that of Mrk 421 (Figure 11), but differed by 
less than 60%. X-ray observations of these blazars showed that their synchrotron peaks are 
vastly different, typically ~ 0.4 to 8 keV for Mrk 421, and well above 100 keV for Mrk 501. 
In light of their vastly different peak energies in the synchrotron region of their spectrum, 
the overall similarity in the value of their peak energies in the Inverse Compton (TeV) region 
is extremely interesting. 

The observed spectrum of Mrk 421 during a period of intense spectral variability is 
monotonically decreasing with energy at all periods. The absorption corrected spectrum 
shows a dramatically different behavior, characterized by the appearance of a peak in the 
spectrum that monotonically shifts to higher energies as the spectrum evolves to higher flux 
states (Figure 12). The data suggests a shift in peak energy between the low and high flaring 
state, regardless of the EBL spectrum used to correct for the 7-ray attenuation. Figure 13 
depict the evolution of E vea ^ with flux level for representative EBL realizations. 

The blazar H1426+428 is the farthest of the blazars detected so far, and its intrinsic 
spectrum is therefore most attenuated. The combined effects of the strong absorption and 
the uncertainties in its observed spectrum yield an absorption-corrected spectrum with a 
wider range of spectral shape than those suggested in previous investigations. Three of the 
eight EBL realizations give rise to an absorption-corrected spectrum that are very similar to 
the two other blazars, and are characterized by the appearance of a peak in the ~ 1 to 5 TeV 
energy region (Figure 14). However, two EBL realizations imply a luminosity peak below 
400 GeV, and two cases suggests a peak energy above 10 TeV. More precise knowledge of 
the spectrum of this blazar, especially extending the observations into the 100 GeV energy 
regime will provide important constraints on the evolution of the EBL with redshift. 

We also considered the uncertainties in the absolute calibration of the energy of ~ TeV 7- 
rays with atmospheric Cherenkov telescopes, and found that they have an important impact 
on the ability to constrain the EBL by the blazar spectrum. The importance of this effect 
is illustrated in Figure 6, and described in §5. Therefore, it will be extremely important for 
the next generation of telescopes to improve these systematic uncertainties. 

The next generation imaging telescopes HESS (Hofmann et al. 2003) MAGIC (Lorenz et 
al. 2003) and VERITAS (Weekes et al.2002) will cover the sub-100 GeV to 10's TeV energy 
regime, covering the transition region between unabsorbed to the strongly absorbed region 
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of the spectrum. At energies below ~ 200 GeV EBL absorption is essentially negligible for 
a redshift of z=0.03, enabling direct measurements of the intrinsic spectra of Mrk 421 and 
Mrk 501 at these energies. In addition for the strongest blazars GLAST will help to bridge 
the 100 GeV regime with measurements at GeV energies providing an important extension 
into region where EBL absorption is negligible. These measurements will provide additional 
constraints on EBL scenarios, since these spectra will have to fit smoothly with the higher 
energy intrinsic spectra of these sources that will be heavily affected by absorption. 

Acknowledgement: We thank Tanja Kneiske for communicating EBL models in digital 
form, and Martin Schroeder for checking the polynomial coefficients in Table 2. WE also 
thank the referee, Felix Aharonian, for his useful comments and for pointing out an error 
in equation (9) in the original version of the manuscript. ED acknowledges NASA's Long 
Term Space Astrophysics (LTSA-2003) Program for support of this work. FK acknowledges 
support by the US Department of Energy. 



REFERENCES 

Aharonian, F.A. et al. 1997, A&A, 327, L5 

Aharonian, F.A. et al. 1999, A&A, 349, 11 

Aharonian, F.A. 2000, New Astron., 5, 377 

Aharonian, F.A. et al. 2000, ApJ, 539, 317 

Aharonian, F.A. et al. 2002a, A&A, 384, L23 

Aharonian, F.A, Timokhin, A. N., & Plyasheshnikov, AS. V. 2002b, A&A, 384, 834 

Aharonian, F.A. et al. 2003, A&A, 403, 523 

Arendt, R. G., Odegard, N., Weiland, J. L., Sodroski, T. J. et al. 1998, ApJ, 508, 74 

Arendt, R. G., & Dwek, E. 2003, ApJ, 585, 305 

Bernstein, R. A., Freedman, W. L., & Madore, B. F. 2002, ApJ, 571, 56 

Biller, S. D., Akerlof, C. W., Buckley, J., Cawley, M. F., Chantell, M., et al. 1995, ApJ, 445, 

227 

Brinkmann W. et al. 2000, astro-ph0011193 

CambresyL., Reach,W.T., Beichman,C.A., & Jarrett,T.H. 2001, ApJ, 555, 563 

Catanese, M. et al. 1997, ApJ, 487, L143 

Catanese, M. et al. 1998, ApJ, 501, 616 



-26- 



Coppi, P. & Aharonian, F.A. 1999, ApJ, 521, L33 

Coppi, P. 2003, Proc. of "The Universe Viewed in Gamma Rays", University of Tokyo Symp., 
eds. R. Enomoto, M. Mori, S. Ynagita (Universal Academy Press, Inc.: Tokyo), p. 

77 

Costamante, L. et al. 2001, A&A, 371, 512 

Costamante, L., Aharonian, F., Ghisellini, G., & Horns, D. 2003, New Astron. Rev. 47, 677 
(astroph/0301211) 

Costamante, L., Aharonian, F., Horns, D., & Ghisellini, G., 2004, New Astron. Rev. 48, 
469 (astroph/0308025) 

Cambresy, L. et al. 2001, ApJ, 555, 563 

Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Daum, A. et al. 1997, Astropart. Phys., 8, 1 

de Jager, O. C, & Stecker, F. W. 2002, ApJ, 566, 738 

Djannati-Ata'i, A., et al. 1999, A&A, 350, 17 

Djannati-Ata'i, A., Khelifi, B., Vorobiov, S. et al. 2002, A&A, 391, L25 

Dwek, E., & Slavin, J. 1994, ApJ, 436, 696 

Dwek, E., & Arendt, R. G. 1998, ApJ, 508, L9 

Dwek, E., Arendt, R. G., Hauser, M. C, Fixsen, D., Kelsall, T., et al. 1998, ApJ, 508, 106 

Dwek, E., & de Jager, O. C. 2001, in "The Extragalactic Infrared Background and its 
Cosmological Implications", IAU Symp. 204, eds. M. Harwit & M. G. Hauser (San 
Francisco: ASP), p. 389 

Elbaz, D., Cezarsky, C. J., Chanial. P., Aussel, H., Franceschini, A., Fadda, D., & Chary, 
R. R. 2002, A&A, 384, 848 

Finkbeiner,D.P, Davis,M., & Schlegel,D.J. 2000, ApJ, 544, 81 

Finley, J. P. et al. 2001, Proceedings of the 27th International Cosmic Ray Conference, ed. 
M. Simon, E. Lorenz, & M. Pohl, (Hamburg, Germany: IUPAP), 2605 

Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, C. L., & Shafer, R. A. 1998, ApJ, 508, 123 

Fossati, G., et. al. 2000, ApJ, 541, 166 

Gaidos, J.A. et al. 1996, Nature, 383, 319 

Gardner, J. P., Brown, T. M., & Ferguson, H. C. 2000, ApJ, 542, L79 

George, I. M., Warwick, R. S., & Bromage, G. E., 1988, MNRAS, 232, 793 



-27- 



Giommi, P., Padovani, P., Perlman, E 2000, MNRAS, 317, 743 

Gorjian, V., & Wright, E. L., & Chary, R. R. 2000, ApJ, 536, 550 

Gould, R. J., & Schreder, G. 1967, Phys. Rev., 155, 1408 

Guy, J., Renault, C, Aharonian, F. A., Rivoal, M., & Tavernert, J-P. 2000, A&A 359, 419 

Hauser, M.G. et al. 1998, ApJ, 508, 25 

Hauser, M.G. & Dwek, E. 2001, ARA&A, 39, 249 

Hillas, A. M., et al. 1998, ApJ, 503, 744 

Hofmann, W., et al., 2003, Proc. of "The Universe Viewed in Gamma Rays", University 
of Tokyo Symposium, eds. R. Enomoto, M. Mori, S. Yanagita (Universal Academy 
Press, Inc.: Tokyo), p. 357 

Horan, D., et al. 2000, Head meeting, No. 23, 05.03 (HoOO) 

Horan, D., et al. 2001a, AIP Conf. Proc. 578, Gamma Ray Astrophysics 2001, ed. S. Ritz, 
N. Gehrels & C.R. Shrader, (Baltimore, Maryland: AIP), 324 

Horan, D., et al. 2001b, Proceedings of the 27th International Cosmic Ray Conference, ed. 
M. Simon, E. Lorenz, & M. Pohl, (Hamburg, Germany: IUPAP), 7, 2622 

Horan, D., et al. 2002, ApJ, 571, 753 

Jauch, J. M., & Rohrlich, F. 1955, "The Theory of Photons and Electrons", (Cambridge, 
MA: Addison- Wesley) 

Kelsall, T., Weiland, J., L., Franz, B. A., Reach, W. T. et al. 1998, ApJ, 508, 44 

Kneiske, Mannheim, K., & Hartmann, D. H. 2002, A&A, 386, 1 

Kneiske, T. M., Bretz, T., Mannheim, K., & Hartmann, D. H. 2004, A&A, 413, 807 

Konopelko, A. et al. 2003, ApJ, 597, 851 

Krawzcynski, H., et al. 2001, ApJ, 559, 187 

Krennrich, F. et al., 1999, ApJ, 511, 149 

Krennrich, F. et al., 2001, ApJ, 560, L45 

Krennrich, F. et al., 2002, ApJ, 575, L9 

Lagache, C, Haffner, L. M., Reynolds, R. J., & Tufte, S. L. 2000, A&A, 354, 247 

Lorenz, E., et al. 2003, Proc. of "The Universe Viewed in Gamma Rays", University of 
Tokyo Symposium, eds. R. Enomoto, M. Mori, S. Yanagita (Kashiwa, Japan), p. 351 

Madau, P., & Pozzetti, L. 2000, MNRAS, 312, L9 

Malkan, M. A., & Stecker, F. W. 2001, ApJ, 555, 641 



-28- 



Mannheim, K., & Biermann, P.L. 1992, A&A, 253, L21 

Mannheim, K. 1993, A&A, 269, 67 

Mannheim, K. 1998, Science, 279, 684 

Mattila, K. 2003, ApJ, 591, 119 

Matsumoto, T. et al. 2000, in ISO Surveys of the Dusty Universe eds., Lemke, Stickel & 
Wilke (Pringer: Berlin), p96 

Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJ, 397, L5 

Marscher, A.P., & Travis, J.P 1996, A&AS, 120, 537 

Metcalfe, L., Kneib, J.-P, McBreen, B., Altieri, B., Biviano, A. et al. 2003, A&A, 407, 791 

Milliard, B., Donas, J., Laget, M., Armand, C, & Vuillemin, A. 1992, A&A, 257, 24 

Miicke, A. & Protheroe, R.J. 2001, Astropart. Phys., 15, 121 

Papovich, C, Dole, H., Egami, E., LeFloc'h, E., Perez-Gonzalez, P. G. et al. 2004, ApJS, 
154, 70 

Partridge, R.B. & Peebles, P.J.E., 1967, ApJ, 148, 377 

Peacock, J. A. 1999, Cosmological Physics (Cambridge: Cambridge University Press) 

Petry, D. et al. 2002, ApJ, 580, 104 

Pian, E. et al. 1998, ApJ, 492, L17 

Primack, J.R., Somerville, R. S., Bullock, J. S., Devriendt, J. E. G. 2001, AIP Conf. Proc. 
558, eds. F. Aharonian & H. Volk (New York: AIP Press), p. 463 

Protheroe, R.J., Bhat, C.L., Fleury, P., Lorenz, E., Teshima, M., & Weekes, T.C. 1998, Proc. 
25th Int. Cosmic- Ray Conf. (Durban), 8, 317 

Protheroe, R.J. & Meyer. H. 2000, Phys. Lett., B493, 1-6 

Puget et al., J.L. 1996, A&A, 308, L5 

Punch, M., et al. 1992, Nature, 358, 477 

Quinn, J., et al. 1996, ApJ, 456, L83 

Renault, C, Barrau, A., Lagache, C, & Puget J-L., 2001, A&A, 371, 771 

Samuelson, F.W. et al. 1998, ApJ, 501, L17 

Soifer, B. T., & Neugebauer, G. 1991, A J, 101, 354 

Stanev, T., & Franceschini, A. 1998, ApJ, 494, L159 

Stecker, F.W., de Jager, O.C., & Salamon, M.H. 1992, ApJ, 390, L49 



-29- 

Tanihata, C. et al. 2004, ApJ, 601, 759 

Vassiliev, V. V. 2000, Astropart. Phys. 12, 217 

Weekes, T. C, et al. 2002, Astropart. Phys., 17, 221 

Wright, E. L., & Reese, E. D. 2000, ApJ, 545, 43 

Wright, E. L. 2001, ApJ, 553, 538 

Wright, E. L. 2004, New Astr. Rev., 48, 465 

Zweerink, J. et al. 1997, ApJ, 490, L170 



This preprint was prepared with the AAS IATgX macros v5.2. 



-30- 



Table 1. EBL Limits and Detection" 



Wavelength 


vl v (nW m 2 sr 1 ) 


Reference 


0.1595 


> 2.9±8-2 


Gardner et al. (2000) 




< 3Q+ 1 - 1 


Gardner ct al. (2000) 


0.2365 


< s.eiL 


Gardner et al. (2000) 


0.3 


12±7 


Bernstein et al. (2000) 


0.36 


2-o±8:2 


Madau & Pozzctti (2000) 


0.45 


4.6±g:S 


Madau & Pozzetti (2000) 


0.555 


17±7 


Bernstein ct al. (2002) 


0.67 


6.715:2 


Madau & Pozzetti (2000) 


0.81 


8-OloJ 


Madau & Pozzetti (2000) 


0.814 


17±7 


Bernstein et al. (200*) 


1.1 


q 7+3.0 


Madau & Pozzctti (2000) 


1.25 


29±16 6 


Wright (2001) 




54±17 c 


Cambresy et al. (2001) 


1.6 


9.01?-? 


Madau & Pozzetti (2000) 


2.2 


7 q+2^0 

'• y -1.2 


Madau & Pozzctti (2000) 




20±6 b 


Wright (2001) 




28±7 C 


Cambresy et al. (2001 


3.5 


12±3 b 


Wright & Reese (2000) 




16±4 C 


Dwek & Arcndt (1998) 


15 


2.4±0.5 


Elbaz et al. (2002) 




2.7±0.62 


Metcalfe et al. (2003) 


24 


9 7+i-i 

z -'-0.7 


Papovich et al. (2004) 


100 


23±6 


Lagache et al. (2000) 


140 


25±7 d 


Hauser et al. (1998) 




15±6 e 


Hauser et al. (1998) 


240 


14±3 d 


Hauser et al. (1998) 




13±2 e 


Hauser et al. (1998) 


240 - 1000 




Fixscn ct al. (1998) 



a See also Hauser & Dwek (2001, Tables 1,3, 4) and Arendt & 
Dwek (2002, Table 6). Errors represent la uncertainties. 

b Derived using the Wright ct al. (2001) zodiacal light model 

c Derived using the Kelsall et al. (1998) zodiacal light model 

c Based on the DIRBE photometric calibration 

f Based on the FIRAS photometric calibration 
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Table 2. Polynomial representations for the different EBL spectral templates" 



Template 6 



polynomial coefficients: <Xj 



0,N 



A(^m)-rangc 



LHH {a} ={0.954270, 0.209657, -2.03417, 2.72242, 3.48905, -9.74560, 

3.14729, 6.70213, -7.21663, 3.01716, -0.593362, 0.0455855 } 
{6} ={0.936188, 0.501445, -1.43078, -0.921636, 1.81965, 0.194080, 
-0.841639, 0.313789, -0.0354290} 

o.5 x {E"«yP°gi (A)]'' + £oMi°gioWP'} 

MHH {1.46525, 0.366351, -2.56385, -0.0579773, 2.17911, -0.533745, 

-0.466031, 0.226609, -0.0276028} 
HHH {1.78163, -0.152917, -4.14983, 4.92401, 3.60609, -11.2143, 

4.74911, 5.42408, -6.51327, 2.78469, -0.552349, 0.0426156 } 
LLH {0.969853, 0.603730, -2.11970, -2.20808, 2.97575, 1.07215, 

-1.92944, 0.670891, -0.0742040} 
MLH {1.47235, 0.175293, -2.55029, 0.545660, 1.67612, -0.954748, 0.139905} 

{1.37639, 0.122876, -0.136912, -0.650521, -5.81963, 5.43507, 

3.92208, -6.80703, 3.32087, -0.711689, 0.0578321} 
HLH {1.58842, 1.31577, -3.67467, -3.24647, 5.03638, 1.01891, 

-3.03305, 1.34105, -0.235389, 0.0143138 } 

{1.89452, -1.81311, 0.965311, 5.89934, -15.5846, 5.20984 

13.0127, -14.4141, 6.18153, -1.23942, 0.0966528} 
LHL {0.925695, 0.381153, -1.25192, -0.160524, 1.11769 -0.429982, 

0.00277214, 0.0114175 } 
MHL {1.43906, 0.257308, -2.22901, 0.432537, 1.41286, -0.787225, 0.113213 } 

HHL {1.73013, 0.226485, -3.33062, 1.07604, 2.24510, -1.71409, 

0.425683, -0.0355017} 

{1.60828, 0.462253, -2.66388, 0.294411, 1.77810, -0.938433, 0.132448} 
LLL {0.980669, 0.585033, -2.41182, -1.90431, 4.02921, -0.175735, 

-2.44258, 1.93244, -0.740240, 0.148669, -0.0123448} 
MLL {1.49096, 0.667023, -2.85743, -1.44668, 2.75698, 0.393004, 

-1.22980, 0.431022, -0.0463593} 
HLL {1.71098, 0.180354, -2.61936, 0.928249, -0.298039, -0.637540, 

2.12254 -1.56680 0.450176 -0.0457011} 
average {1.20000, 0.533686, -1.81116, -1.28135, 1.94804, 0.655142, 

-1.19371, 0.405981, -0.0437552} 



A > 100 


A < 15 


15 < A < 100 


all A 


all A 


all A 


A < 4 


A > 4 


A < 12 


A > 12 


all A 


all A 


A < 2.2 


A > 2.2 


all A 


all A 


all A 


all A 



a Thc EBL intensity in nW m 



is given by: log 10 {yl v ( A) 



Ef=o«il 



logic, (A)^ 



where A ranges from 0.1 to 10 3 (im, unless otherwise noted, and N is the degree of the polynomial. 
b Tcmplates are abbreviated with L, M, and H. See §3.2 for a detailed explanation. 



-32- 



Table 3. Summary of F-Test Results examining the Confidence in the Exponential Rise 

in the Intrinsic Blazar Spectra 



EBL 




Mrk 421 (HEGRA) 






Mrk 501 (HEGRA) 






Template 


Xi/i 


Xv 2 


P{F X , 1, *i) 


Y 2 


X 2 

Au 2 


P(F X , 1, vi) 


EBL 


LHH 


1.877 


1.366 


0.88 


3.18 


1.84 


0.97 


6 


MHH 


2.794 


1.528 


0.97 


5.377 


1.889 


0.99 





HHH 


1.425 


1.249 


0.68 


2.136 


1.755 


0.79 


1 


LLH 


1.249 


1.244 


0.14 


1.847 


1.797 


0.35 


1 


MLH 


1.333 


1.339 


0.0 


1.895 


1.819 


0.43 


1 


HLH 


1.616 


1.096 


0.91 


2.240 


1.726 


0.85 


1 


LHL 


1.455 


1.216 


0.75 


1.912 


1.404 


0.88 


1 


MHL 


1.616 


1.236 


0.84 


2.194 


1.389 


0.95 


()>> 


HHL 


1.760 


1.196 


0.91 


2.34 


1.34 


0.97 


()<> 


LLL 


1.036 


1.111 


0.0 


1.326 


1.446 


0.0 


1 


MLL 


1.226 


1.106 


0.62 


1.450 


1.418 


0.32 


1 


HLL 


1.118 


1.101 


0.27 


1.301 


1.383 


0.0 


1° 



a F x (l, V2) probabilities were calculated to test the confidence that the inclusion of an exponential rise in the intrinsic 
blazar spectrum improves the x 2 of the fit beyond that expected from random fluctuations in the observational data. 
The confidence limit was set at a probability of 95%. Larger probabilities suggest a high degree of confidence that the 
exponential rise in the intrinsic blazar spectrum is real. EBL realizations responsible for such rise in the spectrum of 
either one of the blazars were considered invalid and labeled " 0" . Acceptable EBL are labeled as " 1" . 

b These EBL realizations are still viable when the 7-ray energy scale is shifted down by 15%, reflecting the uncer- 
tainty in the absolute photon energy calibration with atmospheric Cerenkov techniques. 

c The HLL realization can be further excluded from EBL considerations (see table 6) 
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Table 4. Characteristic of the Viable EBL Realizations" 



EBL 




hot 




1 stars 




J- dust 


I dust/ hot 


realization 


(0.1 


- 10 3 /an) 


(0. 


1-5 /im) 


(5- 


- 10 3 /an) 




HHH 




139.8 




91.5 




48.2 


0.35 


LLH 




58.9 




20.0 




38.9 


0.66 


MLH 




93.6 




54.0 




39.6 


0.42 


HLH 




108.5 




70.1 




38.4 


0.35 


LHL 




58.4 




21.1 




37.3 


0.64 


LLL 




48.1 




20.1 




28.0 


0.58 


MLL 




84.5 




56.0 




28.5 


0.34 


HLL 




112.8 




87.3 




25.5 


0.23 



a All EBL intensities are given in units of nW m~ 2 sr _1 . All EBL real- 
izations listed above gave rise to physically "well-behaved" intrinsic blazar 
spectra. The EBL realization HLL can be further excluded because its 
I dust /hat ratio is lower than that in the local universe. See text for further 
details. 



Table 5. Peak energies E pea k(GeV) of absorption-corrected energy spectra of Mrk 421 and 

Mrk 501. 



EBL 


Mrk 


421 


Mrk 


501 


Mrk 501 










template 


Whipple 


HEGRA 


Whipple 


HEGRA 


HHH 


1200 ± 60 


1241 ± 84 


1778 ± 185 


2460 ± 142 


LLH 


450 ± 51 


800 ± 79 


773 ± 159 


1100 ± 133 


MLH 


777 ± 51 


1048 ± 72 


1136 ± 142 


1543 ± 94 


HLH 


1024 ± 50 


1144 ± 65 


1348 ± 107 


1640 ± 75 


LHL 


331±87 


790±107 


1030 ±249 




LLL 


468 ± 52 


818 ± 75 


778 ± 155 


1172 ± 118 


MLL 


837 ± 52 


1043 ± 70 


1190 ± 138 


1546 ± 89 


HLL 


999 ± 47 


1151 ± 62 


1321 ± 123 


1638 ± 70 
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Table 6. Peak energies E pca k(GeV) of absorption-corrected HEGRA energy spectra of 
Mrk 421 and Mrk 501. The observed photon energies were shifted down by 15%. 



EBL Template 



Mrk 421 



Mrk 501 



LHH 


697 ± 81 


1615 ±340 


HHH 


1130±70 


2109± 103 


LLH 


729 ±65 


1055±98 


MLH 


960±60 


1431 ±74 


HLH 


1081±52 


1559±60 


LHL 


712±85 


1969±360 


MHL 


981±77 


2131±152 


HHL 


1156±68 


2082±96 


LLL 


736±63 


1088±88 


MLL 


969±59 


1448±70 


HLL 


1075±55 


1544±57 



Table 7. The Evolution of Peak Energies (GeV) in the Mrk 421 Spectrum with Flaring 

State for Different EBL Realizations 



EBL 






Flaring 


State 






Realization 


1 


2 


3 


4 


5 


6 


HHH 


85±318 


1055±303 


1105±205 


1340±197 


1490±131 


2078±290 


LLH 


31±90 


448±245 


139±123 


581±165 


613±111 


725±258 


MLH 


4±88 


744±250 


608±148 


925±159 


995±104 


1230±213 


HLH 


332±312 


986±238 


978±133 


1158±145 


1231±92 


1457±173 


LHL 




362±321 




539±263 


583±208 


1345±1952 


LLL 


175 ±172 


446 ±244 


143 ±123 


582 ±164 


615 ±110 


730 ±254 


MLL 


47 ±197 


806 ±248 


706 ±145 


982 ±156 


1054 ±101 


1288 ±203 


HLL 


284 ±320 


960 ±237 


942 ±133 


1133 ±145 


1205 ±92 


1430 ±174 



a Energies in units of TeV 
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Fig. 1.— The observed 7-ray spectra of Mrk 421, Mrk 501, and H1426+428 (left panel). 
The right panel depicts the spectral variability of Mrk 421. The lines are analytical fits 
represented by a power law with an exponential cutoff (see eq. 6) to the data. See §2 for 
more details. 
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Fig. 2. — Template spectra representing different realizations of the EBL with select limits 
and detections of the EBL: (1) black triangles-Gardner et al. (2000); (2) white squares- 
Madau & Pozzetti (2000); (3) black diamonds-Bernstein et al. (2002); (4) black circles- 
Wright (2001), Wright & Reese (2000); (5) white diamonds-Cambresy et al. (2001); (6) 
black square-Elbaz et al. (2002); (7) white star-Metcalfe et al. (2003); (8) black star- 
Papovich et al. (2004); (9) white circle-Lagache et al. (2000); (10) black squares-Hauser et 
al. (1998, DIRBE calibration); (11) white triangles-Hauser et al. (1998, FIRAS calibration); 
(12) small black circles: Fixsen et al. (1998). The data and the uncertainties are summarized 
in Table 1. The heavy dashed line going through the 15 /im point represents an "average" 
EBL spectrum, used here only for illustrative purposes. The shaded area in the top figure 
is bounded by the two EBL model spectra presented in de Jager & Stecker (2002), and the 
one in the bottom figure represents the range of EBL intensities presented by Kneiske et al. 
(2004). 
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Fig. 3. — Left panel: The 'average' EBL spectrum as defined in §3 in the text. The shaded bar 
diagram indicates the different wavelength regions depicted in the right panel of the figure. 
Right panel: The 7-ray opacity of a source located at redshift £=0.03 is shown as a thin black 
line. The shaded curves in the figure represent the contributions of the different wavelength 
regions (depicted in the left panel) to the total opacity. A value of H = 70 km s _1 Mpc -1 
was used in the calculations. 
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Fig. 4. — The measured energy spectra of Mrk 421 (Aharonian et al. 2002a) and the 
absorption-corrected spectra using all the twelve EBL realization depicted in Figure 2 are 
plotted versus energy. Also shown are fits to the observations (open circles) given by a power 
law with an exponential cutoff (eq. 6); a parabolic fit (solid curve, eq. 7) to the intrinsic 
spectrum (filled circles); and a parabolic fit with an exponential rise (dashed curve, eq. 8) 
to the intrinsic spectrum. 
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Fig. 5. — Same as Figure 4 for Mrk 501 data (Aharonian et al. 1999) 
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Fig. 6. — The effect of the 15% shift in the energy of the photons on the intrinsic blazar 
spectrum, and on the ability to reject specific EBL realizations. The figure shows that the 
LHH realization that was previously rejected is now viable when the blazar energy is shifted 
down by 15%. More details in §5 of the text. 
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Fig. 7. — The eight EBL realizations that do not give rise to an unphysical intrinsic blazar 
spectrum (see Table 3) are shown in solid lines. Dashed lines represent the three previously 
rejected EBL realization (LHH, MHL, and HHL) that become viable when the 7-ray energies 
are shifted down by 15%. Symbols for the observational limits and detections are the same 
as in Figure 2. 
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Fig. 8. — The 7-ray opacity for the eight EBL spectra depicted in Figure 6. The lightly 
shaded area in the figure is bounded by the opacities corresponding to the two EBL spectra 
adopted by Konopelko et al. (2003), and the darkly shaded one by those adopted by Kneiske 
et al. (2004). 
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Fig. 9. — Top panel: The measured energy spectrum of Mrk 421 in 2000/2001, based on 
the Whipple data (Krennrich et al. 2001) is plotted versus energy. The solid line through 
the observed spectrum (filled squares) is an analytical fit characterized by a power law with 
an exponential cutoff. The other curves depict the absorption corrected spectra for the 
eight viable EBL realizations. Of these the HLL realization (dotted line) can be rejected on 
spectral grounds (see discussion in §6.1). The LHL realization (dashed line) gives rise to an 
unphysical intrinsic spectrum, which is more evident in Figure 8, but cannot be rigorously 
ruled out by the F-test. Bottom panel: the same as the top panel for the HEGRA data 
(Aharonian et al. 2002a). 



-44- 



10 c 



w 

N 

I 

a 

o 

0) 



2 10* 



T3 






10 1 



i 10 3 
w 

CM 

I 

a 

o 

OjO 

Si 

0) 
N 



O 1Q2 






10 1 




H 1 1 — I I I I 



\ 1 1 1 — I I I I 




Mrk 501 (HEGRA) 



observed flux 



10 



E 7 (TeV) 



Fig. 10. — same as Figure 8 for the Mrk 501 data The Whipple data are from (Samuelson 
et al. 1998) and the HEGRA from (Aharonian et al. 1999). 
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Fig. 11. — The peaks of the intrinsic 7-ray spectrum of Mrk 421 derived for the eight viable 
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observations (filled symbols) or the HEGRA observations (open symbols). The dashed line 
represents the points where E pea k(M.rk 421) = E pea h(Mrk 501). The figure shows that peak 
energies for Mrk 501 are systematically higher than those for Mrk 421. 
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Fig. 12. — The observed spectrum of Mrk 421 during various stages of flaring activity 
is plotted versus energy in the top left panel of the figure. The other panels show the 
absorption-corrected spectra for the eight viable EBL scenarios listed in Table 3. The HLL 
realization of the EBL can be rejected because of the low fractional contribution of the dust 
emission to the total EBL intensity (§6.1), and the LHL realization yields an unphysical 
intrinsic blazar spectrum which, however, cannot be rejected by the F-test. 
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Fig. 13. — The evolution of the peak energy of the intrinsic spectrum of Mrk 421 as a 
function of flaring state. All peak energies were normalized to the peak value of the first 
flaring state, except for LHL, which was normalized to the second one. The peak energies 
are listed in Table 6. 
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Fig. 14. — The observed and the absorption-corrected spectra of the blazar H1426+428 are 
plotted versus energy for the eight viable EBL realizations listed in Table 3. The curves in 
the plots represent analytical approximations to the observed and intrinsic source spectra. 
Observed spectra were fitted by a power law (eq. 6 with Eq —> oo), whereas all intrinsic 
spectra were fit with a parabolic function (eq. 7) and shown as solid lines. However, intrinsic 
spectra corrected for the HLH, LHL, MLL, and HLL realizations of the EBL were better 
fitted by an parabolic function with an exponential rise (eq. 8), and are depicted by dashed 
lines in the figure. 



